Observations, theory and implications of 
thermal emission from gamma-ray bursts 



A. Pe'erI and F. Ryde^ 



^ Space Telescope Science Institute, Baltimore, MD 21218, USA; Riccardo 

Giacconi Fellow. 

Department of Physics, Royal Institute of Technology, AlbaNova, 
SE-106 91 Stockholm,Sweden 



2 



Abstract 

Recent analyses show evidence for a thermal emission component 
that accompanies the non-thermal emission during the prompt phase 
of GRBs. First, we show the evidence for the existence of this com- 
ponent; Second, we show that this component is naturally explained 
by considering emission from the photosphere, taking into account 
high latitude emission from optically thick relativistically expanding 
plasma. We show that the thermal flux is expected to decay at late 
times as Fbb ~ t~^, and the observed temperature as T ~ with 
a « 1/2 — 2/3. These theoretical predictions are in very good agree- 
ment with the observations. Finally, we discuss three implications of 
this interpretation: (a) The relation between thermal emission and 
high energy, non-thermal spectra observed by Fermi, (b) We show 
how thermal emission can be used to directly measure the Lorentz fac- 
tor of the flow and the initial radius of the jet. (c) We show how the 
lack of detection of the thermal component can be used to constrain 
the composition of GRB jets. 



1 Introduction 

Despite many efforts, a clear understanding of the physical origin of the 
photons observed during the prompt emission phase in GRBs is still lack- 
ing. The prompt emission spectra are often fitted with a broken power law 
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model (known as the "Band" function, [T]). However, recent Fermi- LAT 
results show that this fit is inadequate in some cases in which high energy 
photons are observed (e.g., [2]). Even more importantly, the "Band" func- 
tion fit does not, by itself, carry any explanation about the physical origin of 
the observed photons. A common interpretation is that the peak observed 
at sub-Mev range is due to synchrotron emission [3]. However, while the 
synchrotron interpretation is consistent with many afterglow observations, 
this interpretation is inconsistent with the majority of the prompt spectra, 
due to steep low energy spectral slopes observed [1]. 

This inconsistency motivated us to search for an alternative explanation. 
Arguably, the most natural ingredient is a thermal component, that should 
exist in the outflow. In principle, thermal photons can originate either by 
the initial explosion, or by any dissipation of the kinetic energy that occurs 
deep enough in the flow, in region where the optical depth is » 1, so that 
the emitted photons thermalize before escaping the plasma. 

In a series of papers O [H [71 [51 O [lOl [TTl [12l |T3j we have extensively 
studied the contribution of a thermal emission component to the overall 
(non-thermal) GRB prompt spectra. Our research is focused on both the 
observational properties [5l [H [HI [13] , theoretical modeling [71 [HI [10] and 
implications of the existence of this component [HI US]- We give here a brief 
summary of our key results. 

2 Observational clues 

Two main difficulties exist in interpreting the observed spectrum. The first 
is that clearly, in addition to the thermal component there is a strong non- 
thermal part. The second is that the properties of the thermal component, 
(temperature and flux) may be time-dependent, hence its signal is smeared 
in a time-integrated analysis, as is frequently done. In order to overcome 
the second problem, one needs to carry a time- dependent analysis. Such an 
analysis was indeed carried by Ryde (OH]). In these works, a hybrid model 
(thermal -|- a single broken power law) was found to adequately describe the 
prompt spectra of 9 burst over the limited BATSE energy band. Moreover, 
the temperature of the thermal component showed a repetitive behaviour: 
a broken power law in time, with power law index T{t) oc after the 

break time at few seconds. 

Repeating a similar analysis on a larger sample of 56 bursts, we found 
(|llj) that the same repetitive behaviour is ubiquitous. Moreover, in this 
work we also considered the evolution of the thermal flux, and found that it 




Figure 1: (Left, center): Histograms of the late time decay indices of the temper- 
ature (left) and thermal flux (center) in a sample of 56 bursts (taken from [H]). 
Right: theoretical results of the flux decay at late times, showing F{t) oc t~'^ (taken 
from [10]) 

too shows a similar behavior: a broken power law, with index F{t) oc at 
late times. The break time, not surprisingly, is the same within the errors 
to the break time found in the temperature behaviour. Histograms of the 
late time power law indices are shown in figure [TJ 

3 Theoretical interpretation 

As thermal photons originate from below the photosphere, one needs to 
study the properties of the photosphere in a relativistically expanding 
plasma. In such a plasma, Lorentz aberration plays a significant role: for 
example, the photospheric radius strongly depends on the angle to the line 
of sight, 6, via rph{9) oc (^^/3 + T^^), where F is the factor Lorentz of the 
outflow [10] . This strong dependence implies that photons emitted from the 
photosphere at high angles are significantly delayed with respect to photons 
emitted on the line of sight: At°^ ~ 30L52r^^^i^ s, where L is the GRB 
luminosity and Qx = Q/W^ in cgs units is used. 

Calculations of the expected decay of the flux and temperature at late 
times are carried under the assumption that the source (the inner engine) 
terminates abruptly at a given time. Photons emitted off axis are delayed 
("high latitude" emission in optically thick expanding plasma), resulting in 
flux decay at late times. Moreover, due to both weaker Doppler boosting 
and energy losses to the expanding plasma at large radii, these photons' 
observed temperature is also lower. By integrating over equal arrival time 
surface in the entire space, it was shown ([TO]) that the flux decays at late 
times as F{t) cx and the temperature as T{t) oc with a ~ 1/2 — 2/3. 



The results of the theoretical calculations of the flux decay (both numerical 
and analytic) are shown in figure [1] (right). 

Clearly, the theoretical predictions are in excellent agreement with the 
observations. While this by itself does not prove that indeed the photons 
that we see are thermal, we find the agreement between theory and observa- 
tions, as well as the repetitive behaviour seen, two very strong, independent 
indications that we indeed are able to properly identify the thermal emission 
component and discriminate it from the non-thermal part. 

4 Implications 

The existence of a thermal emission component, which, as described, is 
a natural outcome of the fireball model, can bring a breakthrough in our 
understanding of the physics of GRB prompt emission. We discuss here 
three important implications of it. 

The relation between thermal and non-thermal emission. As 
GRB prompt spectra is non-thermal, clearly, in addition to any thermal 
component there is a non-thermal part. Thus, according to our picture, the 
observed spectrum is composed of (at least) two separated ingredients, ther- 
mal component originating from the photosphere, and non-thermal compo- 
nent originating from the kinetic energy dissipation above the photosphere. 
The true properties of this non-thermal part can only be deduced after 
subtracting the contribution from the thermal part. This, however, is not 
an easy task: since thermal photons serve as seed photons to Compton 
scattering by energetic electrons (produced by the dissipation mechanism 
above the photosphere), they contribute to the cooling of these electrons. 
Hence, the resulting non-thermal spectrum (from, e.g., synchrotron emission 
or Compton scattering) depends not only on the properties of the accelera- 
tion mechanism (e.g., the power law) but also on the relative contribution 
of the thermal photons. This can lead to a variety of very complex spectra 
(see [ZIIH]), which depend on the various parameters - the photospheric ra- 
dius, the dissipation radius, magnetic field strength, etc. Interestingly, we 
were able to show that for a relatively large parameter space region, a single 
power law may be sufficient to model the non-thermal part over a limited 
energy range. This result is indeed consistent with the single-power law fit- 
ting of the non-thermal part of the spectrum seen in several recent Fermi 
bursts (e.g., GRB090902B, [2j) 

Measuring the parameters of the outflow. One major advantage 
of identifying the thermal emission, is that its radius of origin is known 



- it is the photosphere. The ratio of the thermal flux and temperature, 
TZ = (F/aT^y/'^ (cj is Stefan's constant) must therefore be proportional to 
the photospheric radius (for photons emitted on the line of sight). In fact, 
due to Lorentz aberration, we showed that TZ oc rph{6 = 0)/r. 

In the classical "fireball" model, the photospheric radius depends only 
on two parameters: the (kinetic) luminosity, and the Lorentz factor. Hence, 
for bursts with known redshift, one can use the three measurable quanti- 
ties of emission at the photosphere (thermal flux, temperature and GRB 
distance) to deduce the values of the three unknowns - the luminositj{3, 
the Lorentz factor and the photospheric radius itself. Moreover, in the 
classical fireball model, the dynamics below the photosphere is fully deter- 
mined by the conservation of energy and entropy. One can therefore use 
the values of the luminosity and Lorentz factor at the photospheric radius 
to deduce the radius at which the initial acceleration began. This radius is 
denoted here as tq. Implementing these ideas, we were able to show that for 
GRB970828, at redshift z=0.96, the terminal Lorentz factor is T = 305 ±28, 
and ro = (2.9 it 1.8) x 10^ cm. The statistical errors in the estimate of the 
Lorentz factor, ±10% are by far the smallest from all the methods known 
today. 

Deducing the outflow composition. GRBs show a wide variety of 
properties from burst to burst. While in some bursts, thermal emission 
component is very pronounced (e.g., in GRB090902B the low energy spec- 
trum is so steep and the peak is so narrow that it is very difficult to find 
an alternative explanation to the peak; see [13] ). in other bursts it is much 
less pronounced. As one example, in the bright burst GRB080916C, the 
"Band" model provides a good fit to the spectrum, up to the highest ob- 
served photons energies, at 13.2 GeV [H]. Opacity arguments can easily 
show that these energetic photons cannot originate from the photosphere, 
but from some (significantly) larger radius H, ^ 10^^ cm. The lack of 
the detection of a thermal component as predicted by the baryonic models 
strongly suggests that a significant fraction of the outflow energy is initially 
not in the "fireball" form. Thus, we found [12j the most plausible alternative 
to be Poynting flux entrained with the baryonic matter. The ratio between 
the Poynting and the baryonic flux in this burst is at least ~ (15 — 20). 

^Note that there is an uncertainty in estimating the kinetic luminosity from the flux. 
This uncertainty can be removed once afterglow measurements are available 

^Note that opacity argument by itself does not give the emission radius; in order to 
obtain that, one needs to specify the relation between the radius and the Lorentz factor. 
The results often used in the literature, r — F^cSt rely on assumed knowledge of the 
variability time St, which is highly uncertain. 



5 Summary 



The existence of a thermal emission component, which, as described, is a 
natural outcome of the fireball model, has a potential to bring about a 
breakthrough in our understanding of the physics of GRB prompt emission. 
It is therefore the subject of an extensive on-going research, from all 
aspects. We stress, that the different physical environment prevents using 
the tools developed for studying the afterglow in the study of the prompt 
emission. 
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